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the  potentially  nonlinear  character  of  electron  interactions 
with  EMIC  waves. 


2.  Inhomogeneity  Parameter 

[3]  For  field-aligned  EMIC  waves,  the  cyclotron  reso¬ 
nance  condition  for  electrons  is 
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[i]  The  rapid  loss  of  radiation  belt  electrons  in  the  main 
phase  of  geomagnetic  storms  is  believed  to  be  aided  by 
EMIC  waves,  and  is  usually  analyzed  with  quasi-linear 
theory.  However,  even  moderate  EMIC  wave  intensities 
easily  cause  resonant  electrons  to  respond  nonlinearly,  with 
drastically  different  results.  We  map  out  the  region  of 
nonlinear  behavior  with  a  single  parameter,  and  show  that 
both  the  direction  and  magnitude  of  scattering  can  be 
estimated  by  analytical  expressions.  The  nonlinear 
interactions  typically  lead  to  adveetion  toward  large  pitch 
angles,  rather  than  diffusion  toward  the  loss  eone.  This  is 
expected  to  reduce  the  overall  loss  rate  and  greatly  affect  the 
distribution  of  trapped  electrons.  Citation:  Albert,  J.  M.,  and 
J.  Bortnik  (2(X)9),  Nonlinear  inleraciion  of  radiation  bell  elecirons 
with  electromagnetic  ion  cyclotron  waves,  Geophys.  Res.  Lett.,  36, 

L 1 2 1 1 0,  doi :  1 0. 1 029/2009GL03  8904. 

1.  Introduction 

[2]  Strong  eleetromagnetie  ion  cyclotron  (EMIC)  waves, 
with  typical  intensities  of  ~l -  10  nT,  are  commonly  present 
in  the  outer  radiation  belts  [Meredith  et  ai,  2003;  Fraser  et 
a!.,  2006].  They  are  excited  by  unstable  distributions  of  ring 
current  ions,  and  can  scatter  both  ring  current  and  radiation 
belt  ions  into  the  loss  cone.  They  may  also  strongly  scatter 
radiation  belt  electrons,  especially  near  the  plasmapause  and 
in  high  density  plumes  [Thorne  et  ai,  2005,  2006;  Shprits  et 
a!.,  2008].  Such  scattering  has  mostly  been  modeled  as 
quasi-linear  diffusion  [Lyons,  1974;  Summers  and  Thorne, 

2003;  Gamayunov  and  Khazanov,  2007;  Jordanova  et  a!., 

1996,  2008;  Miyoshi  et  ai,  2008],  which  indicates  that 
EMIC  waves  with  magnetic  field  amplitudes  of  I  nT,  in 
conjunction  with  whistler  mode  hiss  and  chorus,  give  decay 
time  scales  for  I  MeV  electrons  of  less  than  one  day  [Albert, 

2003;  Li  et  ai,  2007].  Similar  calculations  with  stronger 
EMIC  waves  can  reach  the  strong  diffusion  limit,  where 
particles  scatter  into  the  loss  cone  faster  than  they  can  be 
removed  by  the  atmosphere  [Shprits  et  ai,  2009].  However, 
such  strong  amplitudes  raise  doubts  about  the  validity  of  the 
quasi-linear  framework,  and  a  nonlinear  treatment  seems 
called  for,  as  noted  by  Millan  and  Thorne  [2007].  Such 
analysis  for  test  particles  has  previously  been  developed  and 
applied  to  large  amplitude  whistler  mode  chorus  [e.g., 

Nunn,  1974;  Albert,  2002;  Omura  et  ai,  2008;  Bortnik  et 
ai,  2008].  Here,  we  use  a  similar  approach  to  demonstrate 


(I) 


where  =  eB/mc,  m  is  the  electron  rest  mass,  and  e  is  the 
electron  charge  (in  absolute  value).  Resonances  involving 
multiples  of  the  cyclotron  frequency  do  not  play  a  role  for 
field-aligned  waves.  Assuming  the  wave  is  generated  near 
the  equator  and  propagates  away  from  it  [e.g.,  Fraser  et  ai, 
1996],  the  electron  must  be  “co-streaming,"  i.e.,  directed 
toward  high  latitude.  (This  is  opposite  to  the  usual 
requirement  for  resonance  with  ftcld-aligncd  whistler  mode 
waves.) 

[4]  The  dynamics  of  resonant  wave-partiele  interactions 
have  been  treated  by  many  authors  [e.g.,  Nunn,  1974;  Inan 
et  ai,  1978].  Albert  [1993,  2000]  derived  a  “1  ^  dimen¬ 
sional"  Hamiltonian  for  the  normalized  first  adiabatic 
invariant  /  =  uAp  [lmc)^llQ.c’>  the  wave-particle  phase 
and  distance  s  along  a  field  line,  of  the  form 


K(l,i.s)  =  K^{l,.s)yK^{I,s)smi. 


(2) 


where  A^'o  describes  adiabatic  motion  along  a  field  line  and 
A^'i  captures  the  eO'ect  of  the  resonant  wave.  With  field- 
aligned  wave  propagation,  they  take  the  forni 


Ao  —  —  —  4))  +  \l (I  ~  loY  —  I  “ 


2H,.//^’. 


eBy^lmc 

K\  — — ; - tana. 

kc 


(3) 


where  a  is  the  local  particle  pitch  angle  and  4  is  a  constant. 

[5]  For  fixed  5  the  Hamiltonian  is  analogous  to  that  of  a 
plane  pendulum,  whose  phase  portrait  features  two  classes  of 
periodic  trajectories,  with  a  separatrix  between  them.  The 
frequency  of  small  amplitude  oscillations  is  (A^'i  (]^Kq/0!~)^ 
and  the  island  width  is  [4Kd{(^K^JdP)f^.  Albert  [1993] 
considered  the  timescales  for  particle  motion  within  a  frozen 
phase  portrait  and  for  motion  of  the  island  itself,  finding  that 
the  particle  behavior  was  determined  by  their  ratio. 
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R  = 


O-KojOF 

A,{/)-AoA7.«9/)’ 


(4) 


For  \R\  ^  1,  “linear"  motion  is  expected,  with  an 
effectively  random  value  of  the  phase  at  resonance,  leading 
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Figure  1.  Inhomogeneity  parameter  R  versus  equatorial  pitch  angle  and  energy,  for  a  model  of  2  nT  EMIC  waves  in  a  high 
density  plume  at  Z,  =  4.  The  full  equation  (5)  was  used  for  the  left-hand  plot,  and  the  approximate  equation  (7)  was  used  on 
the  right.  The  triangles  indicate  the  partieles  simulated  numerically.  Nonlinear  interactions  are  expected  for  R  <  1. 


to  diffusion  in  /,  while  \R\  1  leads  to  essentially  nonlinear 

motion,  with  phase  bunching  (with  or  without  phase 
trapping),  and  deterministic  changes  in  /.  R  represents  a 
competition  between  the  strength  of  the  wave,  in  Afi,  and  the 
inhomogeneity  of  the  background  magnetic  field  (and 
plasma  properties),  which  enters  through  the  derivative 
with  respect  to  s  (effectively,  latitude  A). 

[6]  Omura  et  al.  [2008]  recently  presented  a  related 
analysis  for  field-aligned  whistler-mode  chorus  waves, 
working  directly  with  the  resonance-averaged  equations  of 
motion,  and  obtained  a  related  nonlinear  oscillator  equation 
with  the  same  properties.  Repeating  their  derivation  for 
arbitrary  (but  field-aligned)  waves  yields  an  “inhomogeneity 
parameter’'  which  coincides  exactly  with  /?,  and  can  be 
written  explicitly  as 


R  =  I  \ul—^ 

Bw  —  1  vx  a;  I  ^2  Up 

7  g;  vj)  I 

2  Qp  ds  j  ’ 


Viil  I 

c  Qj,  Os 


(5) 


evaluated  at  resonance.  Here  ft  is  the  refractive  index,  kc/uj. 

[7]  For  parallel-propagating  EMIC  waves,  //"  is  given  by 
the  Stix  L  parameter  [Stix,  1962],  which  is  conveniently 
written  as 


1 

4 


(6) 


where  A/ is  the  eleetron-to -proton  mass  ratio  m/m^,  /  indexes 
the  ion  species.  Pi  -  ntilmp,  qi  “  and  Z,  is  the  charge  of 
ion  /.  Quasineutrality  requires  '^iZiqi  =  I.  EMIC  waves 
only  propagate  where  (i^  >  0,  which  occurs  in  frequency 
bands  bounded  below  by  cutoffs  and  above  by  resonances  at 
the  ion  cyclotron  frequencies. 

[8]  We  take  the  ions  to  be  H^,  He^,  and  O^,  so  that  p^  = 
U  pHe  =  4,  /3o  =  16,  with  all  =  I.  Typical  storm-time 
values  of  the  ion  concentrations  are  qn  =  0.77,  r///p  =  0.20, 
qo  -  0.03  [Jordanova  et  al.,  2008].  The  value  of  in 
equation  (6)  will  be  dominated  by  the  term  corresponding  to 


the  appropriate  ion  band.  Making  this  approximation,  using 
a  small-latitudc  approximation  of  the  dipole  magnetic  field, 
and  ignoring  density  variation  in  equation  (5)  yield  the 
estimate 


^  ^  — 9A  me  c  p  cos“  a  I 
2L  Rp  me  sin  a  I  — 

From  now  on,  it  will  be  understood  that  “/?”  refers  to  \R\. 

[9]  Figure  1  shows  /?  as  a  function  of  energy  and 
equatorial  pitch  angle  oq  for  field-aligned,  2  nT  EMIC 
waves  in  the  helium  band,  evaluated  in  a  high  density 
plume  at  Z.  =  4,  with  -  15.  The  frequency  uj/Itt  is 

taken  to  be  about  1 .84  Hz,  so  that  a;/Q/  =  0.96  at  the  equator. 
Such  a  large  ratio  is  necessary  to  bring  the  minimum 
resonant  energy  down  to  about  1  MeV  [Summers  and 
Thorne,  2003;  Albert,  2003].  Equations  (5)  and  (7)  give 
very  similar  results  for  R.  Near  the  loss  cone,  the  inter¬ 
actions  are  expected  to  be  linear  (that  is,  diffusive,  as 
discussed  below).  However,  for  a  wide  range  of  pitch  angle 
and  energy  values,  the  interactions  are  expected  to  be 
nonlinear,  so  that  quasi-linear  diffusion  will  not  be  a  valid 
description. 


3.  Equations  of  Motion  and  Analytical  Estimates 

[10]  To  test  these  ideas,  we  follow  particles  numerically, 
using  the  gyrorcsonance-averaged  equations  of  motion 
following  Chang  and  Inan  [1983].  These  equations  take 
the  form 


c/pil  .  p\  OB 

—r  ^ - sin(/>  -  ^ 

dt  c  m'y  2m'yB  as 

^  p^pi^  OB 

dt  c  \k  m^J  ^  2m'^BDs 

d(f)  _  /^kpii  Qp\  ^  eB„.  Zu;  p^  \  cos  4) 

dt  \m'y  1)  c  \k  m'))  p^ 

dt  m7 
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Figure  2.  Trajectories  of  2  MeV  electrons  with  ao  =  15°,  interacting  with  2  nT  EMIC  waves  in  a  high  density  plume.  Even 
for  these  large  amplitude  waves,  the  interaction  is  fairly  linear,  leading  to  nearly  symmetric  scattering  in  oq  and  in  E. 


In  a  dipole  field,  .v  is  easily  transformed  to  latitude  A.  Both 
the  wave  frequency  ^  and  magnetic  field  strength  as 
well  as  the  background  plasma  density,  arc  treated  as 
constant.  The  particles  were  not  allowed  to  reach  the 
location  of  the  crossover  frequency,  where  the  Stix 
parameter  D  =  0.  and  the  wave  can  switch  from  left  hand 
to  right  hand  circular  polarization  [Stix,  1962].  The  wave- 
particle  phase  (j)  differs  from  ^  in  equation  (2)  only  by  sign 
and  phase  conventions. 

[n]  The  Hamiltonian  analysis  of  Albert  [1993]  yields 
analytical  estimates  of  the  net,  nonadiabatic  change  in  I  due 
to  cither  linear  or  nonlinear  resonant  interactions.  These  arc: 


“A/ —  K 


‘A/\7  =  — 


R> 


K, 


d-Ko/dP 


R<  \. 


(9) 


Here  is  related  to  the  value  of  ^  evaluated  right  at 
resonance,  and  is  effectively  random.  Thus  in  the  linear 
regime,  /  suffers  a  random  walk  with  zero  mean,  or  diffuses, 
while  in  the  nonlinear  regime  the  resonant  value  of  ^  is  well 
defined  and  /  experiences  systematic  drift,  or  advection. 
Corresponding  changes  in  energy  and  ao  arc  [Albert,  2000] 


^  sin*  oo  f  SL 

Aao  -  — — - A/,  — ;■ 

(pf  me)  sin  oo  cos  oo 


-A/. 


(10) 


where  is  the  equatorial  value  of  Since  u;  1^^.,  pitch 
angle  changes  will  be  much  larger  than  normalized  energy 
changes. 


4.  Results 

[12]  First,  the  trajectories  of  24  particles  interacting  with 
a  small  amplitude  wave,  -  5  pT,  were  calculated 
according  to  equation  (8).  The  particles  were  started  at  the 
equator  with  E  =  2  MeV,  ao  =  15°,  and  a  uniform  spread  in 
initial  phase,  For  this  arrangement,  R  at  the  expected 
resonance  location  is  nearly  1000.  After  passing  through  the 


resonance,  both  ao  and  E  arc  sharply  scattered.  The  dis¬ 
tributions  of  Aao  and  AE,  though  not  shown,  arc  symmet¬ 
ric  about  zero  and  depend  sinusoidally  on  the  initial  phases, 
which  should  be  similar  to  the  phases  just  before  the 
resonant  interactions.  Furthermore,  the  scattering  ampli¬ 
tudes  agree  quite  well  with  the  prediction  A///„. 

[13]  Figure  2  shows  a  similar  simulation  with  B,,  =  2  nT, 
which  gives  R  ^  2.4.  This  distribution  of  changes  is  also 
fairly  sinusoidal,  though  not  perfectly  so,  and  the  ampli¬ 
tudes  of  the  changes  arc  well  predicted  by  A//„„  shown  as 
the  double-arrowed  line.  The  maximum  change  in  Aao 
about  10°,  which  is  nearly  enough  to  move  the  particles  into 
the  loss  cone  in  one  bounce.  The  pitch  angle  difftision 
coefficient,  ((Aao)^)/r/,  ~  3  x  10^  day  \  is  in  reasonable 
agreement  with  calculations  based  on  broadband  quasi- 
lincar  theory  applied  to  similar  models  of  EMIC  waves  at 
L  -  4.5  [Li  et  ai,  2007;  Albert,  2008].  The  bottom  two 
panels  show  the  expected,  nearly  sinsusoidal  dependence  on 
initial  phase,  as  well  as  mean  changes  (dashed  lines)  and 
mean  changes  plus  or  minus  one  standard  deviation  (dotted 
lines). 

[14]  In  Figure  3,  the  initial  pitch  angle  of  the  2  MeV 
particles  is  changed  to  ao  =  45°,  for  which  R  =  0.54,  a 
moderately  nonlinear  value.  In  fact,  the  distribution  of 
changes  is  very  different  from  above,  with  most  particles 
tracing  very  similar  paths.  This  is  due  to  the  underlying 
“bunching''  of  wave-particle  phases  as  the  particles 
encounter  and  circumnavigate  the  separatrix  between  trapped 
and  untrapped  trajectories  (sec,  c.g.,  Nunn  [1974],  Albert 
[1993,  2000],  Omura  et  ai  [2008]  for  details).  Some  paths 
arc  exceptions,  reflecting  the  fact  that  R  is  not  very  much 
less  than  1.  Nevertheless,  the  predicted  changes  based  on 
AI/s/i,  indicated  by  the  plus  symbols,  give  good  estimates  of 
the  maximum,  and  most  common,  values.  Repeated  inter¬ 
actions  would  cause  the  particle  distributions  to  drift,  not 
diffuse,  in  pitch  angle  and  energy;  the  cumulative  changes 
would  grow  proportionally  to  time,  not  as  Also  note 
that  the  typical  change  in  ao  is  positive,  away  from  the  loss 
cone.  This  is  expected  to  cause  lower  phase  space  density  at 
the  relevant  pitch  angles,  and  lower  overall  loss  rates,  than 
estimates  based  on  pure  diffusion. 
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Figure  3.  Trajectories  of  2  MeV  electrons  with  ao  =  45®,  interacting  with  2  nT  EMIC  waves  in  a  high  density  plume.  For 
these  particles,  the  interaction  is  nonlinear,  leading  to  phase  bunching  with  systematic  increases  in  ao  and  decreases  in  E. 


[is]  None  of  the  particles  in  Figure  3  exhibit  phase 
trapping,  whereby  particles  cross  the  separatrix  and  spend 
many  phase  oscillations  within  it.  This  is  not  surprising 
since  the  particles  are  moving  away  from  the  equator,  as 
required  for  resonance  with  the  EMIC  waves,  and  so 
experience  increasing  magnetic  field  inhomogencity.  Even 
if  /?  <  1  at  resonance,  R  will  be  increasing,  which  is  not 
conducive  to  trapping  [Albert,  2000].  On  the  other  hand, 
particles  may  start  out  inside  the  separatrix.  Figure  4  shows 
2  MeV  particles  starting  with  ao  -  55®  at  the  equator.  The 
linear  resonance  condition,  equation  (1),  is  satisfied  at  A  = 
1.7°,  where  R  ~  0.2.  About  half  the  particles  behave  non- 
lincarly  in  accordance  with  as  shown  by  the  plus 

marks.  However,  almost  half  the  particles  experience  even 
larger,  and  oppositely  signed  changes,  which  develop  over  a 
sustained  interval  of  nonadiabatic  motion,  as  is  characteris¬ 
tic  of  phase  trapping.  Evidently  the  initial  conditions  are 
near  enough  to  resonance  that  some  particles  complete 
several  phase  oscillations  before  becoming  nonresonant. 
Detailed  examination  of  the  trajectories  verifies  that  for 


the  trapped  particles,  0  oscillates  in  a  limited  range  starting 
at  /  =  0.  Physically,  this  would  occur  if  which  has  been 
held  constant,  grew  rapidly  to  “envelop”  regions  of  phase 
space  containing  nearly  resonant  particles. 


5.  Summary 

[i6]  We  have  demonstrated  that  the  cyclotron  resonant 
interaction  of  relativistic  electrons  with  field-aligned  EMIC 
waves  can  be  highly  nonlinear.  Regions  of  expected  non¬ 
linear  behavior  were  mapped  out  with  the  parameter  R, 
using  a  realistic  magnetic  field  amplitude  of  2  nT.  More¬ 
over,  we  have  shown  that  previously  developed  theoretical 
expressions  for  the  expected  pitch  angle  and  energy  changes 
give  good  estimates  of  the  behavior  seen  in  numerical 
simulation  of  the  particle  motion.  Phase  bunching  without 
trapping,  the  more  commonly  expected  version  of  nonlinear 
behavior,  leads  to  rapid  pitch  angle  increases,  away  from  the 
loss  cone,  with  small  changes  in  energy.  In  combination 
with  diflusion,  this  could  lead  to  major,  complex  effects  on 


Bw(pT)»2000.0  E(MeV)-2.0  PA-55.0 


#0  fe 

Figure  4.  Trajectories  of  2  MeV  electrons  with  ao  =  55®,  interacting  with  2  nT  EMIC  waves  in  a  high  density  plume.  For 
these  particles,  the  interaction  is  nonlinear  and  some  particles  are  initially  phase  trapped,  giving  both  systematic  increases 
and  sustained  decreases  in  ao,  and  oppositely  directed  changes  in  E. 
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both  the  loss  rate  and  the  trapped  partiele  distribution,  Sinee 
EMIC  waves  are  believed  to  play  a  major  role  in  the  loss  of 
radiation  belt  eleetrons  during  storms,  sueh  behavior  should 
be  aeeounted  for  in  future  models.  One  approaeh  would  be 
to  use  the  analytieal  expressions  in  a  eombined  diffusion- 
adveetion  formulation. 
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